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Abstract. The suitable basis functions for approximating periodic function are periodic, trigono- 
metric functions. When the function is not periodic, a viable alternative is to consider polynomials 
as basis functions. In this paper we will point out the inadequacy of polynomial approximation and 
suggest a new set of basis functions. The new set does not suffer from the drawbacks of polynomial 
approximation and by using them one can approximate analytic functions with spectral accuracy. 
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1. Introduction . Let us consider a finite sum approximation of a function f(x) 

1 ■ 1 : /(*) « 5>k&(a0, xs[-L,L]. 

fc=0 

If f{x) and its derivatives are continuous and periodic then the suitable basis functions 
are 

(1.2) Mx) = e<t-*)¥ 

and the error goes to zero exponentially fast. If the function is not periodic, a spec- 
tral rate of convergence can be achieved by using Chebyshev polynomials as basis 
functions. Approximating a function by finite sum of Chebyshev polynomials has 
drawbacks (see, for example [TO]). Despite the spectral rate of convergence, Cheby- 
shev polynomials have peculiar characteristics and therefore do not provide an optimal 
set of basis functions. This peculiarity lies in the behavior of their derivatives. While 
Chebyshev polynomials have uniform behavior in the interval, their derivatives are 
non- uniform. In order to clarify this point let us assume that the interval is [—1,1]. 
Chebyshev polynomials in this interval are defined as 

(1.3) Tk(x) — cos(fcarccos(a;)). 
We have 

(1.4) max |T fc (x)| = 1 

x6[-l,l] 



and this maxima is achieved at h — 1 points. Hence, Tk(x) behaves uniformly. Let us 
consider now Upon defining 

(1.5) x = cos(0) 

we get 

(16) dT k _ dT k de _ fcsin(fcfl) 

dx d6 dx sin(0) 
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Hence, in the vicinity of x = we have 

(1.7) maxl^l«fc 

ax 

while at the end points we have 

(1.8) |^(±l)| = fc 2 . 

ax 

Thus, there is no uniformity in the extrema values of the first derivative. This dis- 
crepancy increases for each additional derivative by a power of fc, namely 

(1.9) \ T k \ X a)\ ^ k j l3 

\T^(x b )\ 

where x a and Xb are extrema points in the vicinity of ±1 and respectively. 
Approximating an analytic function by Chebychev expansion is highly efficient when 
the approximated function exhibits a behavior similar to (??). As in the Fourier 
case, the efficiency deteriorates as the behavior of the function approximated moves 
away from this pattern. The obvious alternative is to look for an approximation space 
where the basis functions and all their derivatives exhibit uniform, or 'almost' uniform, 
behavior. The peculiarity of Chebyshev approximation can be demonstrated also by 
considering interpolation. It is well known that polynomial interpolation in equally 
distributed points is not the right approach. The error function has large gradients 
at the boundaries and in some cases, divergence can occur(Runge phenomenon) [4]. 
Insisting on polynomial interpolation, one should resort to non-uniform distribution of 
points. A set of interpolating points which results in exponential rate of convergence 
(for analytic functions) is 

(1.10) Xi = cos — i = 0, • • • , n.15 

\n J 

We have 

(1.11) \Ax max \ = | cos( J - -)| « - 

In n 



(1.12) \Ax min \ = \l-cos(-)\n-(-) 

n 2 \n/ 

The uneven distribution of points does not make sense in the general case. When 
the function approximated does not have large gradients at the boundaries, there 
is no justification for concentrating points there. This inefficient strategy is more 
pronounced when the large gradients of the function are away from the boundaries. 
A way of quantifying the inefficiency of Chebychev interpolation is in the fact that we 
need ir points per wavelength for resolution, compared to two points per wavelength 
which is "a law of nature" as expressed by Nyquist criteria. It seems reasonable to 
look for a set of basis functions where the relevant interpolating points are evenly, or 
almost evenly, distributed. 

In this paper we present a new set of basis functions. This set is composed 
of non-periodic trigonometric functions. Non periodic trigonometric functions were 
first introduced in [7J. The motivation in that case was to overcome the severe 
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stability restriction which results from using Chebyshev polynomials for space dis- 
cretization while solving time dependent pde's. Many researchers reported signifi- 
cant increases in efficiency while using the approximating method described in [7] 
(e.g. pQ, (5], [6], [9], [H]). In the present paper we would like to deal with these basis 
functions just from approximation viewpoint and to emphasize the advantages of us- 
ing these functions as compared to using Chebyshev polynomials. Approximating a 
general function by a linear combination of these basis functions overcomes the draw- 
backs mentioned above and, for analytic functions, the approximation is spectrally 
accurate. The paper is organized as follows. In Section 2 we describe the proposed 
approximation subspace. The basis functions depend on a parameter p, < p < S, 
and this parameter depends on n ( the size of the subspace) . In Section 3 we carry an 
analysis related to accuracy properties and in Section 4, analysis related to resolution 
properties. The paper is concluded in Section 5 in which we present numerical results. 

2. Non-Periodic Trigonometric Subspace. Without loss of generality we 
will consider functions in [—1,1]. Let us define the subspace 

(2.1) S n =span{V>o,---,^n} 
where 

(2.2) tp k (x) = cos{kpx) ft = 0,2,4,..., < p < |zill 
and 

(2.3) tp k (x) = sin(ftjxc) k= 1,3,5,..., 0<p<|./222 

Now, we would like to define an inner-product and, based on this inner product, a 
set of orthogonal functions. For this purpose we will prove first the following Lemma: 
Lemma 1: 

Let T k be Chebyshev polynomial of degree k then, for k even, we have 

(2.4) cos(kpx) = (-l)^T k (sm(px)) 
and for k odd 

(2.5) sin(fcpz) = (-l)^ i T fe (sin(px))Z500 
Proof: 

The proof is by induction on k. It is easily verified for ft = 0, 1 . 
Assume first that k is even. By the recurrence relation of Chebyshev polynomials 
we get 

(2.6) T k+1 (sin (pa;)) = 2 sin (px) T k (sin (px)) - T k ^ (sin (px)) 
and by induction 

(2.7) T k+1 (sin(px)) = (-1)' (2 sin(px) cos(fcpir) + sin((ft - l)px)) . 
Since 



(2.8) 2 sin(pa;) cos(fcpx) = sin ((ft + 1) px) — sin ((ft — 1) px) 
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we get 

(2.9) T fe+1 (sin(px)) = (-1)4 sin((fc + l)px) 

and the proof of the even case is concluded. 

For k odd, using the recurrence relation and induction we have 

(2.10) T k+ i(sm(px)) = (-l) 1 ^ (2 sin(px) sin(fcpx) - cos((fc - l)px)) . 
Since 

(2.11) 2 sm(px) s'm(kpx) — cos((fc — l)px) — cos((fc + l)px) 
then 

(2.12) T k+1 (sm(px)) = (-1) 4 * 1 cos((fc + l)px) 

and the proof of the odd case is concluded. 

The main result of this Lemma - S n is spanned by polynomials in s'm(px). Hence, 
a valid set of basis functions can be 

(2.13) * fc(x)= T fc (^) 0<fc<n. 

We will prove now that these basis functions are orthogonal under suitable inner 

product. 

Lemma 2: 

Upon defining the inner product 

l 

(2.14) <f,g>= J f(x)g(x)w(x)dxl22 

-l 

where 

(2.15) w(x)= , C ° S{PX) 2 ,123 



sin 



(px) y 



the set of functions 



(2.16) $ fe (.x) =T k ( sm (P x A o < k < nl200 

V s m(p) / 

is an orthogonal basis of S n - 
Proof: 

Since $fe(x) are polynomials in sm(px) then € S n and by change of variables 

(2.17) y =^/100 
v ' y sin(p) 

we get 

(2.18) <^),^W>=^/«|* = 0, ^, 



NON-PERIODIC TRIGONOMETRIC APPROXIMATION 



5 



and the proof is concluded. 
Observe that 

(2.19) lim $ k {x) = T k (x). 

p— >o 

Hence , Chebyshev polynomials approximation can be considered as the extremal case 
when we choose p — 0. 

3. Approximating analytic functions by projection on S n . Let f(x) be 
an analytic function defined in [—1, 1]. If 

n 

(3.i) f n (x) = Y, a ^k(x)m 

k=0 

is the projection of f(x) on S n then 



(3.2) a k 



Since 

l 

to qn / * * ^ sin (p) / T k(v) r] a n sin ( p ) 

(3.3 < $ fc , $fe >= / , dy = p k — 

where /3o = 2 and /3k = 1, 1 < fc, we get 

l _ 

/o /i \ 2 f f(y) T k(y) , ,o n 

(3.4) a fc = — / rfy/30 
where 

(3.5) f(y) = f(g(y;p))tii 

and 

(3.6) g(y,p) = - arcsin (y sin(p)) Z90 

P 

is the inverse of (??). Hence, approximating f(x) by the new set of basis functions 
is equivalent to approximating / (g(y,p)) by Chebychev polynomials. As a result, 
it is sufficient to deal with the latter. Observe that now, due to the singularity of 
g(y,p) at y — ±l/sin(p) , the function approximated has large gradient at the 
boundaries and it is justified to use Chebyshev polynomials. The relevant theory 
which discusses polynomial approximation of functions with singularities outside the 
domain of definition here follows [12] . 

Let K be a bounded continuum in C such that K c , the complement K , is 
simply connected in the extended plane and contains the point at infinity. For such 
K there exist a conformal mapping ^(w) which maps the complement of the unit 
disc onto K c [12] . Let O(y) be the inverse of ^>(w) and 

(3.7) B t = {y: \Q(y)\ = t} ( t > 1)Z38 
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denote the level curves in K c then we have the following theorem: 
Theorem 2.1:Suppose t > 1 is the largest number such that F(y) is analytic in- 
side B t . The interpolating polynomials P n (y) with interpolating points y™ that are 
uniformly distributed on K then satisfy 

(3.8) lim max|F(y) - P n (i/)|» = i/390 

n— >oo y£K t 



Since approximating an analytic function by Chebyshev polynomials is equivalent 
to interpolating the function at uniformly distributed points (e.g. Chebyshev points) 
, the asymptotic rate of convergence can be computed by making use of this theorem. 
For K = [—1,1], the relevant conformal mapping is [8] 

(3.9) @(y) = y±^y 2 -l J311 

f(y) is singular at y = ±1/ sin(p) hence , the largest t is 



(3.10) t= * + J' - 1 = i±^M = cot (|) Z313 

sm(p) W sm (p) sm(p) \2J 



and the asymptotic rate of convergence is 
(3.H) ^=tan(f 

Hence, approximating by the new set of basis functions, the asymptotic accuracy is 
ce where 



./312 



(3.12) e =(tan(|)) n 



Z700 



and c is constant which depends on f(y) but does not depend on n or y. 
By choosing 

(3.13) p = 2arctan(£») ,^991 



where e is the machine accuracy, we eliminate the error which results from the singular 
points y = ±1/ sin(p) and get spectral accuracy. Detailed mathematical proof of 
spectral accuracy is given in [3]. 

The choice of p described above is independent of the function we are approxi- 
mating. Obviously, one can do better by choosing the appropriate p for each function. 
This can be achieved numerically by making use of a minimization algorithm which 
finds the parameter p that minimizes the norm of the error vector defined as 

m n 

(3.14) E = J2\Y, ak (P) " f ^ \ 1990 

3 = 1 fe=0 

where Zj, 1 < j < m , are check points randomly distributed in the interval [a, b]. 

Approximating f(x) via interpolation, ak have to satisfy the following n+ 1 
equations 

n 

(3.15) fl fc $ fc(^) = f( x i) < * < nl91 

k=0 
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where {xi}™ = Q is an appropriate set of interpolating points. Due to the equivalence 
mentioned above, a feasible set of interpolating points is 

(3.16) x l ^g{y t ;p) 1 0<i<nl95 
where 

(3.17) yi = cos ^\ .140 

Considering interpolation at (??), it was shown in [7] that the maximal and 
minimal spacing are obtained close to the origin and at the boundaries respectively 
and hence we have 

(3.18) Ax max = - shr 1 (sin(p) sin (-) ) , 141 



(3.19) A.T mi „ = 1 — - sin 1 (sm(p) cos (—") 

v \ \nJ 



.142 



Hence 

2 

(3.20) lim Ax max = lim Ax min = -.149 

p->§ p->-§ n 

Thus, when p is close to ^ the interpolating points are 'almost' equally distributed. 

Let us now show that when p approaches ? ? t ne non uniformity of the derivatives 
described in section I diminishes. Upon using (??) and 

(3.21) cos (9) = 5!^)/ iooo 

sin (p) 

we get 

(3.22) MuW = MuM = 

ax du ax 

where 



'l - [sinfp) cos(6>)] 2 

(3.23) h(9) = -r^l . )' .108 

sin (p) sm(0) 

We have 

(3.24) lim h{9) = ^-.107 

p->f 2 

Consequently, while in the Chebyshev case we have 

dT k _ dT k d9 _ fc sin(fcg) 
^ 5j ~ d9 dx ~ sin(0) ' 

for $fc we have 

(3.26) lim = -fcsin(^). 

p->f ax 2 

Hence, as p is getting closer to -| , the non uniformity of the first derivative 
diminishes. In a similar way one can show uniformity for higher derivatives too. 



s 
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4. On Resolution. In [7] we had conjectured that the function T m (ay) is re- 
solved by k terms where k = [am] ( it was proven later in [5]). The conjecture is 
based on the following analysis. Let 

(4.1) f(x) = s'm(mx) - 1 < x < 1. 

(similar analysis can be done for f(x) = cos(rirx)) . 
By change of variables y — B ^(p) we §et 

(4.2) f(x) = f(y) = sin(marcsin(ysin(|j))) 
where 

(4.3) m=—.mO 

P 

Since (??), choosing r such that m is an odd number we get 
(4-4) f(y) = (-l) m T m (ay) 

where 

(4.5) a = sin(p). 

Since T m (ay) is polynomial of degree m in y, it can be written as 

m 

(4-6) T m (ay) = J2<Tk(v) 

k=0 

while 

2 f 1 T m (ay)T k (y) 



(4.7) a n k ? = / m \» J K ^' dy c = 2, c k = 1 for k > 1. 

Chebyshev polynomials satisfy the recurrence relation 



(4.8) T n+l (x) = 2xT n (x) - T n _x(x).im 
Hence 

(4.9) » = — [2o^m-i(qy)-Jm- 2 (ay)]r fc (y) d fc > , m > 2 . /300 

Since (??) we have 

(4.10) 2yT m -i(y) = T rn (y) + T m _ 2 (y). 
Substituting this relation in (??) we finally get that the coefficients satisfy 

(4.11) a° = 1, al = 0, a\ = a, 

(4.12) a™ ^aa^-a™" 2 , m > 2, 
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(4.13) af = a(c k _ 1 a™- 1 1 + a™- 1 1 )-a™- 2 , 1 < k < m, m>2. 

Carrying out numerical experiments we have observed that , while k < am, a™ 
oscillates, and once k > am, |a™| monotonically decreases. 

Hence, by using {&k(x)}k =0 > tne maximal k is n and therefore, since (??), we 

get 

(4.14) r < r max 
where 

(4.15) r max = " P 



7r sin(p) 



Hence 

71 

(4.16) lim r max = - 

p->f 2 

which is exactly Nyquist criteria. 

5. Numerical Results. In all the examples below we are comparing two algo- 
rithms for approximating analytic functions on a real interval: 

1. Chebyshev 

2. Npt (Non periodic trigonometric). 

The error presented in the tables is defined as 



(5.1) Er 



100 



\ 3=1 



where / is the exact function we are approximating, / is the approximating function 
( either Chebyshev or Npt) and 

(5-2) y. =a+ y-i)*Z£., l<j<100 

are check points, equally distributed in the interval [a, b}. 

For Npt, the tables contain two sets of results. One set, Nptl, is for the case 
where the parameter p is computed according to (??) 



(5.3) p = 2 arctan(e» ), £ = 10 



-14 



and the second, Npt2, is for the case where the parameter p is computed adaptively 
by minimizing (??). The number in the brackets contains the parameter p in each 
case. 

Example 1 

In this example we approximated 

(5.4) f(x) = l —, r, -1< x < 1. 

y ' M ' 2 + cos(lCte)' ~ ~ 

This function behaves uniformly. The results are 
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n 


Er(Chcbyshev) 


Er(Nptl) 


Er(Npt2) 


30 
60 
120 


3.3192e-2 
4.8959e-4 
1.1495e-7 


1.9403c-2 (0.658) 
7.9423e-5 (1.057) 
3.1167c-10 (1.305) 


1.3420c-2 (1.0993) 
2.0303e-5 (1.4457) 
1.9048e-10 (1.4336) 



Example 2 

In this example we approximated the function 



(5.5) f(x) = x 5 cos (50a;) , -1 < x < 1. 

Due to the x 5 term, the function has large gradients close to the boundaries, never- 
theless, the new algorithm outperforms the Chebyshev approximation as can be seen 
by the results presented in the next table 



n 


Er(Chcbyshcv) 


Er(Nptl) 


Er(Npt2) 


40 
50 
60 


6.0165c-2 
3.7617c-3 
1.6279e-4 


3.5717e-2 (0.840) 
5.4146-4 (0.967) 
4.5186e-ll (1.058) 


3.7053c-2 (0.852) 
2.0321e-6 (1.086) 
3.3845e-ll (1.057) 



Example 3 

In this example we approximated the function 

(5.6) f(x) = e- 30x \ -1<x<1 

and the results are 



n 


Er(Chebyshcv) 


Er(Nptl) 


Er(Npt2) 


10 
20 
40 


4.8234e-l 
2.9417e-2 
2.9475e-6 


4.8138e-l (0.0796) 
2.4545e-2 (0.3939) 
7.3752e-8 (0.8402) 


1.3220e-l (1.5708) 
2.0958c-4 (1.5708) 
4.5169e-14 (1.5708) 



In this case, due to the fact that there is a large gradient in the center of the inter- 
val and that the function is almost at the boundaries, the accuracy is significantly 
improved by choosing the optimal parameter (Npt2) which, in this case, is very close 
tof. 



Example 4 

In this example we approximated the function 

(5-7) f(x)= 7=^=^ , -1<x<1 

V 1-1 — ar 

and the results are 



n 


Er(Chebyshcv) 


Er(Nptl) 


Er(Npt2) 


20 
40 
80 


2.5252c-3 
3.7085e-6 
9.7418e-12 


2.8448c-3 (0.3939) 
2.0681e-5 (0.8402) 
3.3488e-8 (1.1783) 


2.5252c-3 (0) 
3.7085e-6 (0) 
9.7418c-12 (0) 
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In this case , the behavior of the function at the boundaries justifies interpolating 
at Chebyshev points. As expected, the minimization process resulted with p = 
which means choosing Chebyshev points. 

Example 5 

In this example we would like to demonstrate the resolution properties of Npt 
compared to Chebyshev. For this purpose we approximated the function 

(5.8) f(x) = sin(1007rx) + cos(IOOttx), -1 < x < 1. 

The results are presented in the next two tables. 



n 


P 


Er(Nptl) 


220 
240 
260 


1.4248 
1.4369 
1.4471 


3.5265e-l 
4.9448e-7 
3.6805e-9 



n 


Er(Chcbyshev) 


320 
340 
360 


2.4532e-l 
2.6849e-4 
1.6117e-8 



Conclusions: We have presented in this paper a new set of basis functions which can 
be used for approximating a general smooth function, defined on a real interval [a, b] . 
The basis functions depend on a parameter p which is a function of the dimension of 
the approximating subspace. When one fixes p to be zero he gets polynomials. Hence, 
polynomials can be considered as a singular member of the family where p is fixed and 
equal to zero. As described in the paper and verified by numerical experiments, the 
parameter p should depend on the dimension of the approximating space and it ap- 
proaches the other extreme value, ^, as the dimension increases. Besides exponential 
accuracy, the approximating function can be computed efficiently using FFT. 
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